Low-energy excess of vibrational states in v-SiC^: the role of transverse dynamics 
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The study of the effects of the density variations on the vibrational dynamics in vitreous silica is 
presented. A detailed analysis of the dynamical structure factor, as well as of the current spectra, 
allows the identification of a flattened transverse branch which is highly sensitive to the density 
variations. The experimental variations on the intensity and position of the Boson Peak (BP) 
in v-Si02 as a function of density are reproduced and interpreted as being due to the shift and 
disappearance of the latter band. The BP itself is found to correspond to the lower energy tail of 
the excess states due to the piling up of vibrational modes at energies corresponding to the flattening 
of the transverse branch. 

PACS numbers: 61.43.Bn, 61.43.Fs, 63.50.+X 



I. INTRODUCTION 



In this paper we address one of the still open questions 
in physics of glasses: the origin of the excess of vibra- 
tional states in the low energy region with respect to the 
usual Debye density of statesi^ This excess, character- 
istic of every glassy system, gives rise to the presence, in 
the region below a few meV, of a more or less pronounced 
peak, observed both in Inelastic X-ray-Scattering (IXS), 
Inelastic Neutron-Scattering (INS), and Raman spectra. 
This peak is universally referred to as Boson Peak (BP). 
Despite its universality, different interpretations and hy- 
potheses have up to now been proposed both on its origin 
and on the nature of the excess vibrational states. On the 
basis of experimental observations and computer simula- 
tions, these modes have been catalogued as spatially lo- 
calized^ spatially delocalized and propagating^^ spa- 
tially delocalized but diffusive in character X Theoretical 
investigations, on the other hand, have interpreted the 
BP in glasses as arising from the lowest-energy van Hove 
singularity of the corresponding crystal)^ or as the pre- 
cursor of a dynamical instability expected in a disordered 
structure^ or as arising from the high frequency dynam- 
ics of a model glass within a mode-coupling-like descrip- 
tion of Gotzeii and of Theenhaus>i£ 

Particularly in the case of v-Si02 at normal density 
the most widely accepted explanation is that the BP 
originates from the piling up of modes near the first 
van-Hove singularity of the transverse acoustic vibra- 
tional branches These modes should mainly involve rel- 
ative rotations of almost rigid Si04 tetrahedra as early 
pointed out by Buchenau and coworkers^ In addition 
to the ambient pressure data, a much useful informa- 
tion is available in the case of v-Si02 for higher densi- 
ties, obtained both from in-situ measurements and from 
permanently densified samples. In general, densifica- 
tion results in a shift of the BP towards higher energies 
and in a simultaneous decrease of its intensityiiii^i^i^ 



In addition to the experimental data, the BP behavior 
with density has been successfully reproduced by a series 
of simulationsiiiiiLiS Recently, very accurate INS mea- 
surements on densified vitreous silica (d-SiO^) have con- 
firmed these effects iS2i2ii2S In densified samples the BP 
lies at higher energies than in vitreous silica at normal 
density and has a rather weaker intensity. In densified 
samples the inelastic signal in the BP energy range shows 
a more marked dependence on the scattering wavevector 
Q as regard its intensity and peak position. We present 
here simulation results on v-Si02 at different densities 
which, when compared with the existing experimental 
data, can help in clarifying the origin and the nature of 
the excess modes as well as the intensity and shift effects 
on the BP as a function of density. 



II. SIMULATION 

The simulated system was a glass consisting of 680 
Si02 molecules interacting through the modified BKS 
two-body interaction potential) 2 ^* with periodic bound- 
ary conditions. After a long thermalization at a den- 
sity of 2.2 g cm -3 and at 5000 K the system was slowly 
cooled down to 300 K at constant volume. At this tem- 
perature, an energy minimization was performed to lo- 
cate the equilibrium configuration at room pressure. The 
box size was then stepwise scaled by 1.5% amounts up 
to the final density of 4.0 g cm -3 . Once reached this 
value, the density was decreased in the same way. At 
each step the energies of the individual realizations were 
minimized using geometrical methods to obtain the rel- 
ative equilibrium configurations. This procedure mimics 
a real hydrostatic pressure experiment at room tempera- 
ture as described in detail elsewhere^ The pressure was 
estimated by using the virial theorem and compared to 
experimental data>^2^ After the density cycle the zero 
pressure realization had a density of about 2.8 g cm" 3 . 
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The dynamical properties of the two zero pressure config- 
urations, and of the higher density ones, were computed 
in the harmonic approximation by diagonalizing the dy- 
namical matrix. In particular, we have computed the 
density of states g(E) and the longitudinal and transverse 
currents spectra (C(Q,E)) which, in the one-excitation 
approximation, are given by: 

c t (Q> E ) = -jms; ' • g f W) ex p o • R ^ 2 

■ S(E - E p ) 

Cl {Q,E) = -*&L= ■ S p |S„(Q x e p (n))exp(z Q ■ R n )\ 2 

■S(E~E p ) (1) 

where a, (3 indicate Si and O, and Q—Q/\Q\- 

The current spectra in transverse polarization are, con- 
trary to the longitudinal currents, not experimentally 
measurable, and thus they are not directly observable 
by means of inelastic scattering experiments at least in 
the low-Q and \ow-E ranges. They are nonetheless very 
useful in understanding the details of the vibrational dy- 
namics, as will be explained in the following section. 



III. DISCUSSION 
A. Current Spectra 

Usually the dynamical data obtained via simulation 
are reported showing the dynamical structure factor, 
S(Q,E). This is mainly because this procedure has the 
advantage of being directly comparable to the experimen- 
tal spectra obtained, for instance, from standard INS and 
IXS measurements. On the other hand, the current spec- 
tra, defined as C(Q, E) = E 2 S(Q, E)/Q 2 , are much more 
useful when one is interested in characterizing the vibra- 
tional spectra. In fact, each current spectrum, taken at 
a fixed E and C(Q,E), gives the Q-cnaracterization of 
the individual vibration with energy E, being the cur- 
rent strictly related (in harmonic approximation) to the 
power spectrum of the Fourier transform of the corre- 
sponding eigenvector. The computed current spectra, 
both for longitudinal and transverse dynamics, are shown 
in Fig. ^ upper and lower panel, respectively) as a con- 
tour map plot for v-SiC"2 at p = 2.2, 2.8, 4.0 g cm" 3 . 
The most evident feature is that in the longitudinal cur- 
rents a periodicity is present, even if much less evident 
than in a crystal. The same feature is present in the 
transverse ones, but with a longer quasi period which 
extends beyond the limit of the figure. 

As discussed elsewhere jliLSSi this difference is due to 
the fact that the longitudinal dynamics is affected by the 
average height of the SiCU tetrahedron, which is the main 
structural unit of v-SiC>2, while the transverse vibrations 
are sensible to the average height of SiCU tetrahedron 
decorated with oxygen atoms. In the low-Q longitudinal 




FIG. 1: Computed currents relative to neu- 

tron scattering grey scale maps for v-Si02 at 
three different densities (from left to right): 
p = 2.2 g cm" 3 (P = GPa), p = 2.8 g cm" 3 (P = GPa 
densified); p — 4.0 g cm" 3 (P = 35 GPa). Upper panels: 
longitudinal currents; lower: transverse. The current spectra 
have been divided by Q 2 for visualization purposes. 



spectrum, in addition to the quasi-periodic pattern, two 
well observable features are present in the two lowest- 
density samples, which consist of protunding peninsu- 
lae, nearly parallel to the Q axis, centered at about 20 
and 100 meV. These features have been interpreted^ as 
a trace of the vibrational dynamics of the crystobalite 
structure (the crystalline counterpart of v-SiOg), which 
shows in that energy range flat dispersion bands. It 
should be noted that the plots of Fig. ^ are shown as 
a function of the modulus of the wavevector, which is a 
relevant quantity in glasses, thus loosing every informa- 
tion on the Q-direction. If this convention were used also 
in representing the dispersion curves of a crystal, differ- 
ent Q-directions in the Brillouin zone would superimpose 
giving qualitatively the same effect as in Fig. In the 
Q-E range between 5 and 15 nm _1 and 10 to 25 meV 
(where the flat band is observed in the longitudinal cur- 
rents) a much more evident flattening is present in the 
transverse currents. This coincidence has suggested the 
interpretation that the flattened band in the longitudi- 
nal currents arises from a spilling of transverse dynamics 
into the longitudinal one, due to disorder. 

The maxima of the longitudinal and of the transverse 
currents in the constant Q cuts, are reported in Fig. [21 
for three different densities. 

In the starting configuration, p = 2.2 g cm -3 , open 
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FIG. 2: Main maxima of the longitudinal (open symbols) 
and transverse (crossed symbols) current spectra, for three 
different densities. Squares: p = 2.2 g cm -3 (P — GPa); 
diamonds: p = 4.0 g cm~ 3 ; circles p = 2.8 g cm~ 3 (P = 
GPa densified). 



FIG. 3: Vibrational density of states of v-Si02 at different 
densities: full squares: p = 2.2 g cm -3 ; full circles: p = 
2.8 g cm -3 ; up triangles: p = 3.2 g cm -3 ; down triangles: p = 
4.0 g cm -3 ; open circles: compressed sample p — 2.8 g cm -3 
(zero pressure after pressure cycle). 



and crossed squares of Fig. the main maxima relative 
to the longitudinal currents have a rather defined disper- 
sive character up to very high energies (about 100 meV), 
while the transverse dispersion one flattens at small Q- 
values. The transverse maxima remain centered, within 
the uncertainty, at the same energy where a minor peak is 
observed also in the longitudinal current spectra, whose 
intensity is roughly 50 to 80% of the main featured By 
increasing the density up to 4.0 g cm -3 , the local struc- 
ture of the glass, which at normal density consists of 
fourfold coordinated Si ions, is almost completely sub- 
stituted by an octahedral one (typical of stishovite, the 
crystalline form of silica, stable at high densities). As a 
consequence, also the dynamical properties change and 
this can be seen in Fig. [21 Both longitudinal and trans- 
verse maxima follow a nearly linear dispersion law up to 
the highest energies and no flattening is observed. 



B. Boson Peak 

From our simulations we calculated also the density of 
states as a function of the density. The result are shown 
in Fig. for selected densities. The low energy part of 
the total density of states shifts toward higher frequencies 
and decreases in intensity with increasing density. As a 
consequence, the BP strongly decreases in intensity and 
shifts at high energies^. 

The reduction of the intensity of the BP with the den- 
sity can be ascribed to two concurring effects. The first 
one is an overall reduction of excess states, induced by 
the pressure, and the second one is the shift towards 



high energies of the excess which, being weighed by an 
l/E 2 factor, appears as having a lower intensity. Actu- 
ally, both mechanisms are effective in v-SiC>2 as we will 
discuss later. In Fig. 0] we present the maxima of the 
transverse currents computed at several densities. In the 
p = 2.2 g cm~ 3 sample, after a linear increase a flatten- 
ing of the maxima of the transverse currents is observed 
for q- values higher than 6 nm -1 , at about 15 meV. This 
plateau progressively shifts, and eventually disappears 
with increasing density. In the highest density sample it 
is no more observed and only a linear dispersion relation- 
ship is present. By decreasing the density the plateau is 
at least in part recovered, as shown in the upper panel 
of Fig. 0] It is worth noting that in the densified sample 
at P = GPa (p = 2.8 g cm -3 ) the flattening is present, 
but it saturates at an energy of about 5 meV higher than 
in the starting realization, indicating an hysteretic be- 
haviour. 

The shift and the disapperance of the plateau help in 
interpreting the intensity and the energy effects on the 
BP. The excess of states shifts to higher energies and 
progressively disappears. The fact that these effects take 
place in an energy region, i.e. 15 meV, higher than the 
one where the Boson Peak is observed, about 4 meV, is 
only an apparent contradiction. As a matter of fact, the 
presence of the l/E 2 factor severely affects the spectral 
shape of the excess states by weighting only their low- 
energy tail. Indeed, the true excess of states should be 
observed only in the difference between the computed 
density of states and its Debye approximation. 

The determination of the actual shape of the vibra- 
tional excess is not a simple task, since it implies the 
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FIG. 4: Lower panel. Maximum of the transverse current 
spectra computed at different densities with increasing den- 
sity. Diamonds, circles, triangles and squares correspond to 
p — 2.2 g cm -3 , p = 2.8 g cm -3 , p = 3.2 g cm -3 , and 
p = 4.0 g cm -3 , respectively. Upper panel. Maximum of 
the transverse current spectra measured at different densities 
with decreasing density. Down triangles, circles, up triangles 
and squares correspond to p — 2.4 g cm -3 , p — 2.8 g cm -3 , 
p = 3.2 g cm -3 , and p = 4.0 g cm -3 , respectively. The 
crossed symbols refer to the zero pressure realization, before 
and after the pressure cycle. 



knowledge not only of the density of states, but also of 
the relative weight of the crystalline one, or at least its 
Debye approximation. The density of states can be ob- 
tained from inelastic neutron scattering data, and also, 
to a certain extent, from unpolarized Raman data. The 
Debye density of states can be independently estimated 
via Brillouin light scattering, while heat capacity exper- 
iments give the proper normalization factorsi^ 

The result of this procedure is shown in Fig. [5] where 
INS-deduced g(E) is shown with the relative estimate of 
the Debye contribution. The excess of states (full line 
in the figure), obtained by subtracting the Debye from 
the DOS, is a very broad unstructured band centered at 
energies higher than 7 meV. By dividing this latter by 



FIG. 5: (a) Dotted line: INS density of states of v-Si0 2 ; 
dash-dotted-line: Debye approximation; full line: vibrational 
excess, i.e. density of states minus Debye; dashed line: den- 
sity of states divided by E 2 . (b) Same as (a) but for the 
unpolarized Raman spectrum of v-SiC>2. Here the full line 
represents the density of states weighted by the Raman cou- 
pling function. 



E 2 its maximum shifts towards lower energies, less than 
4 meV, as is shown in Fig. [S] dashed line. This latter 
peak corresponds to the Boson Peak observed both in 
INS and IXS measurements of the dynamical structure 
factor. Actually, the observed BP lies on the low energy 
tail of the true excess of states, since the contribution 
of the remaining high-energy vibrations is almost sup- 
pressed by the 1/E 2 factor. A similar effect is observed 
also in the depolarized Raman spectra. Here one does 
not observe directly the density of states, but this quan- 
tity weighed by the Raman coupling function^ as shown 
in Fig. EJb). A qualitatively similar behaviour is found 
even if the scattered intensity is somewhat deformed by 
the Raman coupling function C(E). 

The difference between the position of the maxima in 
the excess of states as found in simulations (about 17 
meV as in Fig.^J) and in experiments (about 7 meV as in 
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Fig-El, can derive from two main causes. From the com- 
putational point of view, the discrepancy may derive from 
the potential used. Moreover, the experimental density 
of states is obtained through some approximations, like 
f.i. the incoherent approximation and the estimate of the 
multiphonon contribution. Also, in the simulations the 
excess states could be overestimated due to unphysical 
quenching times currently used in Molecular Dynamics, 
as pointed out recently by AngelL 3 - Finally, the behav- 
ior of the transverse current maxima with density shown 
in Fig. 01 not only accounts for the intensity decrease 
of the BP, but also explains the coherence effects of the 
BP measured with in INS and IXS as a function of the 
density. The BP in normal v-SiC>2 has a small coher- 
ence effect, i.e. it is rather insensitive to the exchanged 
wavevector^ On the contrary, densified samples have a 
marked coherent behavior at least at small Q-values^S 
From Fig. 0] full circles, one sees that the coherence ef- 
fects in the BP range, i.e. variation in position and/or 
intensity, are expected in the linear dispersing range in 
the case of the p = 2.2 g cm -3 sample, being the disper- 
sion curve flat for wavevectors higher than 6 mn" 1 . In 
the densified sample at zero pressure, p = 2.8 g cm -3 , the 
linear region extends up to 9 mn" 1 . In the latter sample 
it is indeed possible to study, using IXS, a more extended 
linear region and the direct observation of a Brillouin 
peak directly in the energy cuts is possible in spite of the 
resolution limitations. Hence, in the densified samples 
our simulations predict (i) in the low-Q, low- £7 range a 
dispersing region where a Brillouin peak is observed ex- 
tending to wavevectors higher than than in normal silica 
measurements, (ii) the appearance of the BP (nearly Q 
insensitive) at wavevectors (and corresponding energies) 
higher than in v-SiC>2, (iii) a BP intensity weaker than 
that of the normal density sample. All these predictions 



are confirmed by the experiments of Foret and cowork- 
ers^ while the energy shift and the intensity effect have 
also been reported by Inamura and coworkers^ 



IV. CONCLUSIONS 

In conclusion, the comparison between experiment and 
simulation allows us to affirm that the BP in v-SiC>2 cor- 
responds to the low-energy tail of a flattened transverse 
branch, whose center lies at relatively high energies. This 
is also in agreement with earlier studies ji^i which as- 
cribed these vibrations to hindederd librations of coupled 
SiC>4 tethrahedra. The effect responsible for the occur- 
rence of the BP which we have described here is valid 
for v-SiC>2, and it cannot be excluded that other mech- 
anisms are active for different systems, like the idealized 
ones studied in Ref. Nevertheless, it suggests a gen- 
eral approach for the study of the BP in glasses. It is well 
known that every kind of disorder induces a sort of level 
repelling and a consequent piling up of modes at the ex- 
tremes of the vibrational spectrum, as shown by Elliott 
and coworkers, (see for instance Ref. 32 and references 
therein) . 
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